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Abstract 

The small x behaviour of the non-singlet structure function is studied within the double 
logarithmic approximation (DLA) of perturbative QCD. Since there is neither k t nor 
0 ordering in the ladder Feynman graphs, the predicted non-singlet quark densities for 
the HERA kinematical range (x ~ 10~ 3 ) exceed the values calculated from the small-x 
approximation of the conventional Altarclli-Parisi evolution by a factor up to ten. 
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1 INTRODUCTION 


The small x behaviour of the non-singlet structure functions in deep inelastic scattering 
plays an important role for the precise description of the quark densities. Indeed, in order 
to check the Gottfried sum rule or to estimate the fraction of the initial nucleon spin 
carried by the quark (the spin crisis problem) one has to extrapolate the data into the 
small x region, and therefore one has to know the behaviour of the structure functions 
at small x. 

The kinematical region x -C 1 is of the Regge type, hence there are two possibilities 
for its investigation: 

(i) to use conventional phenomenological formulae of the Regge type for the structure 
function description at the hadron level (see e.g. the review Q and references therein); 

(ii) within the framework of the perturbative QCD, to calculate and sum up the 
Feynman graphs yielding the leading contributions to the structure functions at the 
partonic level. 

In the latter case one should take into account the Feynman graphs to all orders in 
the coupling constant for x -C 1. The most popular way to carry this out consists of 
extrapolating (see e.g. i, i- i) the Altarclli-Parisi equations, which describe DGLAP 
evolution of the structure functions with respect to ln(Q 1 2 ) DU, H, 0 into the region 
x <C 1. This equation sums up in the region x ~ G(l) all the leading contributions 
coming from the Feynman graphs of the ladder type for the kinematical region where 
the transverse momenta kjT of virtual partons (quarks and gluons) are strongly ordered 

hr < h T < ... < k nT < t/q 2 . ( 1 ) 

The quantity Q 2 = — q 2 > 0 denotes the square of the photon mass. 

Indeed, it has been proven in [0 that an integration over kj in the region (1) yields 
the main contribution to the singlet and non-singlet structure functions at x rs./ 0(1). 

However, such ordering is not valid for summing up the leading contributions when 
iCl. In this case the well known result [] is that the condition (1) should be replaced 
by the ordering of the longitudinal Sudakov variables cq, A for momenta A as follows: 


Oil -c «2 -C ... C «n < 1, 

1 > A > ... » /3 n . 


( 2 ) 


Integration over /q in the region (2) yields the leading contributions to the structure 
functions, in particular, the double-logarithmic (DL) contributions for the non-singlet 
structure functions (NSSF). Thus, one can evaluate NSSF by the direct calculation of the 
Feynman graphs of the ladder type. However there is another, simpler way to do this, 
namely, to make use of the Infrared Evolution Equation (IREE) for NSSF, i.e. of the 
evolution equation with respect to the infrared cut-off in the transverse momentum space. 
This method has first been applied in |10] to the description of the elastic scattering of 
quarks, and it has then been generalized to the investigation of inelastic scattering in 


1 Firstly, this result has been obtained in 

QED. See also the review [pi. 


for the investigation of the process e + e 


+p~ 
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QED and QCD (see [1C], [111 ,[12] and references therein).On the whole, this method 


has turned out to be an efficient way to calculate scattering amplitudes at high energies. 

In the present paper we calculate the non-singlet structure functions in the double 
logarithmic approximation (DLA) by constructing and solving both the integral equa¬ 
tion of the Bethe-Salpeter type and IREE for it. In contrast to calculations based on 
extrapolating the conventional DGLAP into the region j; C 1, we take into account both 
the contributions ~ [a s ln(x ) ln(Q 2 )] n and the DL-contributions ~ ( a s ln 2 x) n in every 
order in a s . The latter ones become even more important at very small x. 

The paper is organized as follows. In Section 2 we review the main features of 
calculating Feynman graphs in DLA. We derive the linear integral equation for NSSF, 
taking into account the running QCD coupling constant a s (Q 2 ) within DLA. In Section 
3 we construct and solve IREE for NSSF. Here we assume a s to be fixed. The solution 
is given in terms of the Mellin transformation. In Section 4 we discuss our results and 
compare them with the conventional expression for NSSF obtained with Altarelli-Parisi 
evolution equation. 


2 INTEGRAL EQUATION FOR THE NON¬ 
SINGLET STRUCTURE FUNCTION 

In this section we give the outline of the Feynman graph calculation in DLA developed 
firstly in [|J ( see also the review |J) for the forward e + e _ —> p + p~ scattering in QED 
and use it for constructing the integral equation for the non-singlet structure function 
fNs{x, Q 2 ) (Jns = q n .s.( x i Q 2 )) i n deep inelastic lepton-nucleon scattering (DIS) at small 
x. In the framework of the perturbative QCD f^s is given by the graph in Fig. 1 where 
the dashed lines correspond to the virtual photon with momentum q and the straight 
lines correspond to the quarks. The wavy lines we assign for gluons. The blob in Fig. 
1 means that the radiative corrections in DLA are taken into account. The dotted 
line crossing the blob means that the intermediate s-channel particles are on-shell. We 
investigate Jns in the kinematical region 

|i 2 <<<3 2 « 2 pq (3) 

i.e. at small x, where x = Q 2 /2pq, p denotes the initial quark momentum and p is the 
infrared cut-off. The leading (the double-logarithmic) contributions to f^s come from 
the ladder Feynman graphs (Fig. 2). Non-ladder graphs, though yield DL-contributions, 
do not contribute to f^s because they compensate each other totally. The reason for 
such a compensation in QCD and QED is that all non-ladder gluons are nearly on-shell 
and can be treated like the bremsstrahlung in the forward scattering, cancelling each 
other due to the destructive interference. 

Thus, DL-contributions to the non-singlet structure function come from the ladder 
Feynman graphs in Fig. 2 only. To evaluate them, one should solve the linear integral 
equation in Fig. 3. for the sum of DL-corrections, A, to the elastic forward quark- 
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antiquark scattering amplitude A(t,/3) (see Fig. 3): 


A(t, (3) 



(4) 


where a s is the QCD-coupling constant, C F = ( N 2 —1)/2 N for the colour group SU(N). 
We have used in (4) the Sudakov parameterization for k in the following form 

k = aq' + (3p + kp, —kp = t > 0, (5) 


where 




q 2 Pv 

2 pq 


( 6 ) 


and the analogous parameterization for k' . The nonhomogeneous term, Aq, in (4) is 


A 0 (t,/3) =4nC F a a (t/P). (7) 

The limits in the integral over 0 in (4) corresponds to the /5-ordering (see e.g. 0); 

0> (3. ( 8 ) 

The lower limit in the integral over t is given by some infrared cut-off /i, originated by 
applicability of the perturbative QCD to the problem , while the upper one comes from 
the condition k 2 = a(3s <| t | (where s = 2 pq). Indeed, if k T k T the value of a in 
(5) is a ~ ( k T — k' T ) 2 /0s ~ t' / s0 and, when t' t0 / f3, the longitudinal part of k ij 
becomes so large, 

= sa(3 ~ t (3/(3 > t, (9) 

that it destroys the logarithmic structure of the integrand dt'/t' in (4) (see |5] for the 
details). 

So, within the DLA we put 

t m = t0/(3 . (10) 

Let us notice that for small x (i.e. when 0 (3) this value is larger than the upper limit 

when ti are ordered. Now let us show that the argument of the QCD coupling constant 
in (4) is also equal to t m . To find it, one has to consider the quark loop inserted into the 
gluon propagators in Fig. 3 (see also Fig. 4). The integral over the invariant mass of the 
loop looks as / dM 2 /M 2 . This logarithmical behaviour continues up to M 2 

rS ' / ^771 = t(3 03 

where the large value of M 2 corresponds again to the essentially large k jj ~ M 2 (3j(3 ~ t. 
Thus, the argument of a s should be t m Q. 

Finally, to obtain the non-singlet structure function, f N s, we integrate A over t: 


Ins{x,Q 2 ) = e 2 q \5(l-x) + JL/ A{t,0) j } , 


( 11 ) 


2 Note that for small x (i.e. when [3 (3) this value is much greater than the argument in the 

canonical Altarelli-Parisi equation (t m = t). 
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with s m = s/4 and 


(3 — x + t/s . 


( 12 ) 


e q in (11) denotes the electric charge of the initial qnark. Here the upper limit, s/4, 
is fixed formally by the kinematics, but, in reality, the bulk of integration in (11) is 
determined by dependence on f3. At 1 > Q 2 , the term t/s increases the value of (3 
crucially, and as the amplitude A(t, /3) falls down with an increase of (3 the essential 
value of t in equation (11) is of the order of Q 2 . It should be noticed that the expression 
(11) corresponds to the integration with the transversely polarized photon i.e. to cr T . 
In the case of the longitudinal heavy photons one has to write, in the denominator, Q 2 
instead of t 


L f A(t,f3)dt/t 
cr _ jP_ 

° T ]'A(t,l3)dt/Q 2 ’ 

id 


(13) 


where f3 in (13) is given by (12). 


3 IREE FOR THE NON-SINGLET STRUCTURE 
FUNCTION 


In this section we construct and solve the evolution equation for NSSF with respect 
to the infrared cut-off in the transverse momentum space. It is more convenient to 
consider, instead of the non-singlet structure function, f^s, the scattering amplitude, 
M, for the forward quark-photon scattering where photon is off-shell (see Fig.4). The 
Optical theorem gives the simple relation between f^s and M. 

Indeed, 


47 t a 


2 tr{pe(p + q)e} „ 1 


( 2 >„r hs = A mM 


(14) 


where e /t is the photon polarization vector and a is the QED coupling constant. As we 
are discussing the forward scattering, the amplitude M depends on s, Q 2 and on the 
masses of quarks taken into account. 

There are graphs with the ultraviolet and the infrared singularities among the Feyn¬ 
man graphs contributing to M. We can easily avoid the problem of the ultraviolet 
singularities by noticing that we will treat any single Feynman graph contribution to M 
as the result obtained with the dispersion relations from its imaginary part (with respect 
to s) which has not the ultraviolet singularities. In order to get rid of the infrared sin¬ 
gularity problem, one should introduce some infrared regularization. Let us define the 
infrared cut-off parameter, p: 


/i <C k lT , i = 1,2,... (15) 

where Ayr denotes the momentum of virtual particle i (a quark or gluon), transverse to 
the plane formed by momenta p and q' (see (5)): 

kirP = k iT q = 0. (16) 
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Such a regularization does not destroy the gauge invariance, so we can fix any gauge for 
the virtual gluons at convenience. In the present paper we use the Feynman gauge. 

It was shown in |TT|] that /i, being the minimal transverse momentum, also takes part 
of a new mass scale. Thus, we can neglect masses of the virtual quarks and to be free of 
the infrared singularities in the same time. In DLA, the value of /i can be chosen much 
greater than A qcd- 

As M now depends on s, Q 2 , and on the new mass scale /i, in the form 


M = M(s/y 2 ,Q 2 /y 2 ), 


(17) 


one can differentiate it with respect to fi. Let us notice that 



dM dM 
dz dy 


(18) 


where 

z = ln(s/n 2 ), y = ln(Q 2 /y 2 ). (19) 

Expression (18) is the left-hand side of the IREE for M (the left-hand side of the equation 
in Fig.5). 

Now let us obtain the right-hand side of the IREE. In the first place, let us notice 
that the ordering conditions (8) do not mean the ordering for fj in the ladder Feynman 
graphs. Indeed, integration of such graphs over cq yields (see [9;]) 


< t 2 /P 2 < • • • < t n /(3 n . 


( 20 ) 


Since U are not ordered, integration over any of t t gets y 2 as the lowest limit. However, 
the region (20) can be considered as a superposition of the regions where the transverse 
momenta are ordered, with f,; [i — 1, • • • ,n) being the minimal transverse momenta in 
the region fi*. Therefore, in every such region fb, y is the lowest limit only for integration 
over kiT whereas integration over other transverse momenta does not involve y at all. 
In other words, the ladder Feynman graph ( and therefore the sum of them) can be 
treated as a convolution of two ladder graphs, so that, in each the graph, integration 
over transverse momenta begins from ti in the region i.e. for them k % x takes place 
of /i, being a new infrared cut-off when the quark masses are neglected. Let us denote 
such a minimal transverse momentum as kj*. Ordering (20) means that the quark pair 
with the minimal kj* can be in any place of the ladder graph in Fig. 2. In the case of 
tj-ordering the quarks with the minimal k t would be in the fixed place of the ladder 
graph (in the bottom of the graph in Fig. 2). 

Thus, beyond the Born approximation, the amplitude M can be represented as the 
convolution of two amplitudes (the last term in the right-hand side of the equation in 
Fig. 5). The first of these amplitudes is the same amplitude M with replacement y by 
kx and s by 2 kq. The second one is the amplitude Mq for the forward quark-antiquark 
scattering. 

Adding the Born amplitude Mb to the right-hand side of the equation, we obtain the 
equation in Fig. 5 for f^s i n the graphic form. The letters inside the blobs in Fig. 5 
denotes the infrared cut-offs. This equation can be easily rewritten in the analytic form 
in the momentum space with using the standard Feynman rules. 
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It is convenient however, to perform the Sommerfeld-Watson transformation which 
partly coincides with the Mcllin transformation : 


A+ioo 


M(z,y)= I ^(s/y 2 )“F(u,y) 


( 21 ) 


A—200 


where A is chosen so that the integtation contour in (21) was on the right of the sin¬ 
gularities of F. The variable iv corresponds to the angular momentum in the complex 
momentum plane. Having performed the Mellin transformation, we obtain the equation 
for F in the following simple form: [] 


d ly = lO"' F 


ojF + = zr-z.fn(u)F (22) 

where /o is the Mellin amplitude for the forward scattering of on-shell quarks in DLA 


101. There is not a contribution coming from the Born amplitude Mb in (22) because 


M b does not depend on // in the region (2). 

As equation (22) is linear, one can easily obtain its solution: 

F = #(w) exp{(-u) + /o/8ir)y} 


(23) 


<h(u;) in (23) should be specified with a certain boundary condition. Let us fix it at the 
point y = 0: 

F(uj,y) 1^0= (24) 

Thus, in order to specify $ we should construct and solve IREE for the amplitude of 
quark-photon forward scattering, M, where the photon is (nearly) on-shell. Immediately 
we obtain for its Mcllin amplitude, F(u), similar equation: 


coF = C 0 + —f^)F(u) 

_ -47 rtr{pe(p + q)e}e 2 q 

On — - 

2 pq 


(25) 


(26) 


The first term in the right-hand side of (25) is the contribution of the Born graph which 
depends on y when y — 0. 

Solution of (25) gives us 

F = Co/(u> — /o/87t 2 ). (27) 

The Mellin amplitude fo(u>), firstly obtained in |10;], satisfies the same differential 
equation like F{u) with obvious replacements F by f 0 in equation (25) and C 0 by a 0 , 


a o — Q 2 Cf 

where g is the QCD-coupling constant. Thus, the equation 

w/o = a 0 + / 0 2 /(8ir 2 ) 


(28) 


(29) 


3 For details of the method see 


and the references therein. 
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for /o is non-linear, though very simple. Its solution is 


fo = An 2 uj[l - sj 1 - a 0 /(27r 2 o; 2 )] (30) 

where the sign minus before the square root sign was chosen so that 

fo U>i~ do/uj (31) 

where ao/cu is the Born approximation for fo. Eq. (29) permits to simplify (27) because 
it gives 

u - /o/(8vr 2 ) = ao/fo (32) 

Thus, combining (29), (27) and (23), we obtain 

F(v,y) = ^r fo(uj) exp{-[u - f 0 /(8TT 2 )}y}. (33) 

It leads to 

M( f^ ) = 7ki ^(^ToMexpWoMAS^)}. (34) 

X—ioo 

With using (14) and the fact that in DLA 

ImM tv —ndM / 'dln(s) (35) 

we obtain for the non-singlet structure function the following expression 

2 A+ioo 

= fk / (36) 

X—ioo 


4 DISCUSSION 


The expression (36) for the non-singlet structure function takes into account DL con¬ 
tributions to all orders in the coupling constant, and therefore can be represented as 
an expansion in powers of a s . In terms of the Mellin variable u, we can consider the 
expansion of fo into a power series in 1/u and integrate the result over u according to 


— (-) u u~ k - 1 
2ni v x ' 


- k k k 


k = 1, 2,---. 
As 


fo ~ 47T 2 (^- + t) + 0(u~ 7 ) 


'2ca 8w 3 16cn 5 ' 

with a = 2a s Cp/7T, we obtain from (36) 


fNS 


4vr 2 e 2 
g 2 Cp J 2m W 

exp 


du /1 \ u / a, 

- 1 -) K— 


+ 


a 


+ 


GT 


{?/[ 


2u 8cu 3 16u 5 
2 


a a 
+ 


+ 


a 3 , 


4 u 16cn 3 32cn 5 1 


(37) 


(38) 


(39) 


7 



It should be stressed that for the non-singlet function at small x one deals with the 
expansion (a s /to 2 ) k , which is much more singular ( at to —> 0) than the expansion 
( a s /to) k typical for the singlet (BFKL-pomeron) case. Thus, the summation of the 
whole DL series becomes very important. 

If to a, one can neglect all terms in (39) proportional to ( a/to 2 ) k (k > 1), compared 
to the first one. It leads to the well known expression /vs for the non-singlet structure 
function: 


\+ioo 


Ins ~ Ins = el 


A—200 


did / 1 \ c 
2iri '-x' 


exp( 


Auo' 


(40) 


Expression (40) is obtained from the conventional Altarclli-Parisi evolution equation 
which is valid in the region x ~ 1 where it correctly sums up all the leading InQ 2 contri¬ 
butions to /vs, by retaining in the anomalous dimension only the singular part in 1/c o. 
However, as it is clear from (39), it fails in the region x -C 1, where DL-contribntions, 
proportional to ln 2k [l/x) becomes extremely essential and therefore should not be ne¬ 
glected. 

To demonstrate this let us expand (36) and (40) in inverse powers of to and integrate 
them over to with using (37). We obtain 

Ay- = 5(x - 1) + y[/n( l/x) + ln(Q 2 /p 2 )\ + 


Ins 


+—[2ln 3 (l/x) + 6ln 2 (l/x)ln(Q 2 //i 2 ) + 
9b 

+3ln(l/x)ln 2 (Q 2 / /i 2 )] + 0(a 3 ), 
2 

= S(x - 1) + jln(Q 2 //i 2 ) + ^/n(l/ x)ln 2 (Q 2 / /j 2 ) + 0(a 3 ). 


(41) 

(42) 


Though the first, the Born terms, in (41), (42) coincide, higher loop corrections are quite 
different. Indeed, for the ratio of the first loop contributions (terms proportional to a) 
we obtain 

_ ln(l/x) /t . 

ri = + /n(g 2 //i 2 ) ’ ( } 

and for the ratio of the second loop contributions (terms proportional to a 2 ) 


r 2 = 1 + 


2ln(l/x) 2ln 2 (l/x) 


+ 


ln(Q 2 /n 2 ) 3ln 2 (Q 2 / /i 2 )' 

If we put, for example, 

yu 2 = 0.1 GeV 2 , Q 2 = 20 GeV 2 and x = 10 -3 , 


(44) 


(45) 


we obtain 

ri 2.3, r 2 ~ 4.7. (46) 

Thus, being equal in the Born approximation, /vs and the conventional non-singlet 
structure function /vs differ a lot when the radiative corrections are taken into account. 
This difference increases with the order of perturbation expansion. 



Now let us estimate the ratio fws/f ns at small x, using their asymptotic expressions. 
With the standard technique (the saddle point method), we find that at x -C 1 

~ fNs ~ ^0 1,A ^ h ^ p) (47) 

where we have used the notations p = — ln(x ); b = \/a, y = In (Q 2 /p 2 ). 

The asymptotic expression for f^s at the same region is 

fNS ~ l/^ 1 “ ^)\Z5 exp{6(p + ?//2_?/2/4p)} ' (48) 

Thus, at x —> 0 

Ins/Ins ~ ^(d/f/) 3/4 exp{6( v /py -p-y /2 + ?/ 2 /V)} (49) 

For example, if we choose the same values (45) for /r 2 , s, Q 2 and put aq = 0.3, we obtain 
from (49) the estimate /ns/Ins ~ 7.89. More accurately, a numerical evaluation of fj^s 
and /ns given by (36) and (40) yields the following result: 

W Ins ~ 5.56 . (50) 


Our results for other values of x and p, with Q 2 = 20 GeV 2 , are presented in Table 
1. The value of a s for p 2 = 0.1 was put by hands equal to 0.3 whereas for other /j 
a s = a s (p 2 ). 


Table 1. 

/i 2 , GeV 2 

X 

R = f ns t Ins 

a s 

0.1 

0.001 

5.56 

0.3* 

0.1 

0.01 

2.94 

0.3* 

0.1 

0.1 

1.64 

0.3* 

1.0 

0.001 

10.26 

0.394 

1.0 

0.01 

4.52 

0.394 

1.0 

0.1 

2.10 

0.394 

4.0 

0.001 

12.71 

0.283 

4.0 

0.01 

5.94 

0.283 

4.0 

0.1 

2.75 

0.283 


In the Mellin representation, the Altarelli-Parisi equation for f^s can be written as 

^Ins 


d ln(Q 2 


— —cn + g C F /ijj} f N s 


(51) 


IREE gives the Q 2 -dependence of /ns as: 


9fNS 

dln(Q 2 


= -u + fo(u)]f NS . 


(52) 
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The first terms in the brackets in (51), (52) are not essential. They disappear if one uses 
the Mcllin transformation in the form 



instead of ( 21 ). 

Thus, the only essential difference between (51) and (52) is the difference in the 
’’kernels” in these equations. When x ~ 1, the replacement of / 0 by its Born value 
makes (52) coincide with (51). In a sense , the opposite replacement g 2 Cp/u; by fo(u>) 
corresponds to taking into account all next-to-leading corrections of the type of (a s /u ; 2 ) n 
to the Altarclli-Parizi kernel. Being small in the region u ~ 1, they become large when 
u —> 0 (x —> 0). In this connection, let us notice that these corrections are proportional 
to ( a s /u) 2 ) k only due to the fact that integration over all transverse momenta yields large 
contributions from the region of small hr (h t > /i). Hence, if one wants to take into 
account non-perturbative contributions to f^s ( see e.g. 0 ), one has to consider them, 
besides the leading order, in the non-leading orders also. 

Thus, in the present paper we have shown that the applicability of the conventional 
expression (40) for the non-singlet structure function at the small x region is quite 
doubtful. The standard and widespread Monte-Carlo programs (e.g. HERWIG and 
LEPTO) should be changed when they operate with the /ns in the region x <C 1. 
The main reason for the descrepancy is that in the Regge kinematical region x <C 1, 
the transverse momenta of the virtual quarks in the ladder graphs are not ordered, in 
contrast to the transverse momentum ordering in the Altarclli-Parisi evolution at x ~ 1. 
It leads to replacement of expression (40) by /ns given by (36). 

Such an amendment leads to a drastic change of the asymptotic behaviour of the non¬ 
singlet structure function. Indeed, the conventional formula (40) leads to the behaviour 
for x —> 0 

Ins ~ exp J 2as ° F ln(Q 2 //i 2 )ln-, (54) 

V 7T X 

while the IR evolution equation gives at x —> 0 





(55) 


thus demonstrating the more singular behaviour of /ns compared to f^s at small x. 
Such a behaviour can be written as 


Ins ~ x 7 . 

where 7 is an effective power. One can see that, for f N s, its value 


7 = 


2 (x s Cf ln(Q 2 /jj 2 


\ 7 T ln(l/x) 
is smaller at very small x than 7 from (55) which predicts 

7 = \/ 2 a,C F / 7 r. 


(56) 


(57) 


(58) 


10 












ACKNOWLEGEMENT 


We are grateful to J.Bartels and to L.N.Lipatov for useful discussions. We would like 
to gratefully acknowledge the hospitality of the DESY Theory Group and the financial 
support of the Volkswagen Stiftung. 


11 



References 

[1] M. G. Ryskin. Surveys in High Energy Physics 7 (1994) 137. 

[2] G. Ingelman, ’’Physics at HERA” edited by W. Buchmiiller, v.3, p.1366; Hamburg, 
1991. 

[3] G. Marchesini, B. R. Webber, Nucl. Phys. B310 (1988) 461; 

B. R. Webber, ’’Physics at HERA” edited by W. Buchmiiller, 

G. Ingelman, v.3, p. 1354; Hamburg, 1991. 

[4] S. Forte and R, D. Ball. CERN-TH.7436/94 

[5] V. N. Gribov and L. N. Lipatov, Sov.Journ.Nucl.Phys. 15 (1972)438 and 675. 

[ 6 ] G. Altarclli and G. Parish Nucl.Phys. 126 (1977)297. 

[7] Yu. L. Dokshitser, Sov.Phys.JETP 46 (1977)641. 

[81 V. G. Gorshkov,V. N. Gribov, L. N. Lipatov and G. V. Frolov. Sov.J.Nucl.Phys. 

6 (1967) 95. 

[9] V.G.Gorshkov, Sov.Phys. Uspekhi 16 (1973) 322. 

[10] R. Kirschner and L. N. Lipatov, Nucl.Phys. B213 (1983) 122. 

[11] B. I. Ermolaev and L. N. Lipatov. Int. J.Mod.Phys. A4 (1989) 3147. 

[12] P. Chiappetta and B. I. Ermolaev. Int. J.Mod.Phys. A 8 (1993) 5353. 


12 



Figure captions 

Fig.l. Graph for fjyg- The blob indicates that the radiative corrections in DLA are 
taken into account. 

Fig.2. Graph of the ladder type for f^s • Crosses on the lines mark on-shell intermediate 
particles. 

Fig.3. Equation of the Bethe-Salpeter type for A(t,/3). 

Fig.4. The amplitude M ( // denotes the infrared cut-off). 

Fig.5 The evolution equation for M and M. 
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